AN  IBM-65()  PROGRAM  FOR  THE  COMPUTATION 
OF  THERMAL  GRADIENTS  IN  MASS  CONCRETE 
STRUCTURES  , 

f  , 

E.  E.  McCoy,  et  al  ’ 

I  _  ' 

Army  Engineer  Waterways  Experiment  Station 
Vicksburg,  iv!ississippi  > 

April  1960 


National  Technical  Informatieii  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road.  Springfield  Va.  22i51 


TECHNICAL  REPORT  NO.  6-540 
April  I960 


AN  IBM -650  PROGRAM  FOR  THE 
COMPUTATION  OF  THERMAL 
GRADIENTS  IN  MASS 
CONCRETE  STRUCTURES 


TECHNIC-ftL  REPORT  NO.  6-540 
April  I960 


U.  S.  Army  Engineer  Waterways  Experiment  Station 
CORPS  OF  ENGINEERS 

Vicksburg,  Mississippi 


\/ICKSBURG  MISS 


PREFACE 


In  Cjvil  VJorks  engineering  des_gn  studies  cons.'def’ation  is  ai'' ays 
given  to  the  possible,  undesirable  effects  of  the  therrcal  gradieni  -  which 
may  develop  in  mass-concrete  structures  during  construction.  Par.  t  ’aph 
5-O3  of  EM  lllC-2-2000,  Standard  Practice  for  Concrete,  includes  reference 
to  two  analytical  methods,  the  Carlson  method  and  the  Schmidt  m.ethod_  which 
may  be  used  to  compute  temperature  changes  in  concrete. 

The  investigation  leading  to  the  preparation  of  the  comput"*r  program 
described  in  this  report  was  authorized  by  letter  fi'om  the  Office,  Chief  of 
Engineers,  dated  21  March  19^8,  subject  "FY  195^  Civil  Works  Investiga¬ 
tions  -  Assignment  of  CW  Item  9^6."  The  work  was  accomplished  as  a  joint 
project  of  the  Concrete  Division,  under  the  supervis‘’on  of  Messrs.  Thomas  B. 
Kennedy,  Bryant  Mather,  and  E.  E.  McCoy,  and  the  Technical  Services  Divi¬ 
sion,  under  the  supervision  of  Messrs.  C.  B.  Patterson  and  D.  L,  Neumann. 
This  report  was  prepared  by  Messrs.  McCoy  and  Neumann.  Pvt.  Stephen  W. 
Closs,  Jr.,  assisted  during  the  developmental  stages  of  this  work  and 
provided  se’''eral  routines  which  are  included  in  the  final  program. 

Col.  Ecunund  H.  lang,  CE.  was  Director  of  the  Waterways  Experiment  Sta¬ 
tion  during  this  investigation  and  preparation  of  this  report.  Mr.  J.  B. 
Tiffany  was  Technical  Director. 


Preceding  page  lilanit 


1 1 1 


CONTENTS 


Papjje 

PREFACE .  iii 

SLWARY .  vii 

PART  I:  INTRODUCTION .  1 

The  Problem .  1 

Theory  and  Method . 1 

Thermal  Properties  of  Concrete  and  Other  Required  Data  ...  2 

PART  II:  COMPUTER  PROGRAM .  5 

The  Problem  and  Method  of  Solution  .  5 

Description  of  Equipment  .  5 

Input-Output  Data  Format  .  5 

Operatin(f  Instructions . 8 

Flow  Die.gram .  10 

PART  III:  EXAMPLES .  ih 

Check  Problem .  14 

A  Construction  Problem .  l8 

Accuracy .  21 

Concluding  Remarks  .  ...  22 

REFERENCES  . .  24 


Preceding  page  blank 


V 


SUMMARY 

A  computer  program  for  the  computation  oi  the  unsteady  state  temper¬ 
atures  that  develop  during  the  construction  of  a,  concrete  dam  in  layers 
(lifts)  is  described.  Information  pertainii'g  tc  the  following  ambient  and 
structural  conditions  and  thermal  properties  is  necessary  for  a  computa- 
^  tion;  ambient  air  temperature  and  foundation  tempeiature,  concrete  temper¬ 
ature  at  time  of  placement,  lift  nhicKness,  number  of  days  between  lifts, 
and  adiabatic  temperature  rise  and  thermal  diffusivlty  of  the  concrete. 

The  computation  may  be  continued  for  fifty  lifts. 

The  computation  is  based  on  a  nianerical  procedure  for  the  solution 
of  the  heat  flow  equation  for  one  dimension.  The  space  ■'.'.terval  is  one- 
quarter  thickness  oi  a  lift,  and  the  time  interval  is  one -quarter  day. 

The  program  is  vrritten  for  the  basic  IBM-65C  computer  wita  \  2CCC-word 
dnrni. 


Examples  are  gi'^en  to  illustrate  the  types  of  problems  for  which  the 
program  is  suitable.  In  one  example  the  maximum  temperature  for  ^-day 
placement  of  lifts  was  6  degrc*es  less  tha.n  that  for  3-day  placement,  and 
the  maximum  temperature  for  7-day  placement  was  5  degrees  less  than  tnat 
for  5"day  placement.  In  another  example  the  advantage  of  use  of  re¬ 
frigerated  concrete  Is  illustr.ited.  The  maximum  tempei’ature  computed  in 
this  example  was  only  h  degrees  higher  than  the  initiaJ  foundation 
temperature . 


Preceding  page  blank 


AH  IBM- 6 50  PROGRAM  FOR  THE  COMPUTATION  OF  THEBMkL  GRADIENTS 
IN  MASS  CONCRETE  STRUCTURES 


PART  I:  INTRODUCTION 

The  Problem 


1.  Extensive  studies  of  the  effects  of  temperature  changes  in  a  dam 

have  seldom  been  undertaken  because  to  compute  the  temperatures  manually 

entails  so  much  labor,  and  the  chance  of  error  is  so  great,  that  a  thorough 

study  by  manual  computation  is  impracticable.  Also,  since  comparative  data 

are  desired  from  such  studies,  it  is  advantageous  to  investigate  as  many 

conditions  as  may  have  an  important  effect  on  the  temperatures  in  a  dam, 

but  the  work  is  miiltiplied  by  each  set  of  conditions  that  is  investigated. 

6* 

Although  the  problem  of  adverse  thermal  behavior  has  increased  over  the 
past  tw'o  or  three  decades,  principally  because  of  increased  sizes  of  dams, 
a  thorough  study  of  all  uhe  important  placement  conditions  which  could  af¬ 
fect  the  temperatures  in  a  dam  has  never  been  made  for  even  a  single  proj¬ 
ect.  In  contrast,  on  completion  of  development  of  the  computer  program 
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described  herein  it  was  possible  to  undertake  studies  for  a  dam  involving 
l6  different  combinations  of .concrete-placement  scheduling,  climatological 
factors,  placement  temperature  of  the  concrete,  and  heat  liberation  of  dif¬ 
ferent  types  of  cement,  at  a  cost  estimated  to  be  less  than  one-tenth  that 
of  a  manual  computation. 


Theory  and  Method 

2.  The  computation  of  lu.steady  state  temperatures  in  a  dam  is  most 
conveniently  resolved  into  the  problem  of  computing  temoeratures  at  a  suf¬ 
ficiently  large  number  of  poiiits  and  at  su-Pficiently  rapid  intervals  to 
yield  the  desired  accuracy.  Any  method  of  computing  the  temperatui’e  at  a 
point  in  a  solid  material  or,  better,  the  average  tempei’ature  of  a  small 
unit  vol'me  inclosing  the  point,  must  take  into  consideration  inflow  of 


*  Raised  numbers  refer  to  References  at  the  end  of  this  report. 
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heat  to  the  point,  outflow  of  heat  from  the  point,  and  the  net  generation 

of  heat  due  to  chemical  or  physical  processes  active  in  the  material  at  the 

point.  Two  methods  have  been  extensively  used  during  the  past  three 

decades  for  the  computation  of  temneratiires  in  a  structure.  They  are  the 
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Schmidt  metho  I.,  which  has  been  -widely  applied  to  problems  invol-vLng  -vari¬ 
ous  structures  and  shapes,  and  the  Carlson^  method,  which  was  developed 
specifically  for  concrete  walls  and  dams.  The  latter  method  ^vas  used  for 

the  program  described  in  this  report.  This  method  is  "essentially  a  step- 

^  2 

by- step  integration,"  or  convenient  "step  method"  for  the  numerical 
solution  of  the  heat-flow  equation  for  one  dimension  with  added  term  for 
heat  generation  within  the  concrete.  Analyses  in  one,  two,  or  three 
dimensions  can  be  made,-**-  although  the  procedure  for  only  one  dimension  is 
given  in  the  paper  by  Carlson. 

3.  The  temperatures  of  primary  interest  in  a  dam  are  in  the  lower 
lifts.  Because  of  the  thiclmess  of  the  d.'am  at  the  elevations  of  these 
lifts,  temperatures  in  them  can  be  effectively  controlled  by  variation  of 
construction  practices.  Since  the  thick  lifts  are  relatively  self- 
insulated,  lateral  heat  flow  may  be  considered  negligible,  and  a  one- 
dimension  analysis  is,  therefore,  sufficient.  The  Carlson  method  for  one 
dimension  (the  vertical),  exactly  as  described  in  reference  1,  was  used  in 
the  development  of  the  computer  program. 

Therms-.l  Properties  of  Concrete  and  Other  Required  Data 

4.  Computer  personnel  must  be  furnished  the  necessary  Information 
per-taining  to  each  computation.  This  information  is  noimaliy  specified  by 
or  supplied  by  an  engineer  who  is  familiar  with  construction  practices  and 
with  the  determination  of  thermal  properties  of  concretp*.  The  following 
conditions  define  a  problem: 

a.  Ambient  temperat’ora .  A  different  value  of  ambient  air  tem¬ 
perature  Ta_  may  be  given  for  each  period  of  time  beiween 
lifts  (see  d  on  the  following  page)  but  not  for  a  shorter 
period  of  time.  A  constent  value  may  be  used  throughout. 


**  See  -x’eferenee  j,  a’athor's  c-Iosui’e. 
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b.  Concrete  temperatuje.  The  temperature  of  the  concrete  T^ 
at  the  time  of  placement  must  be  gi/en.  It  may  be  any  value 
either  colder  or  varmer  or  equal  to  the  ambient  air 
temperature . 

c.  Lift  thickness ♦  The  program  can  handle  any  reasonable  lift 
thickness  (4X);  however,  the  thickness  must  be  held  con¬ 
stant  throughout  a  computation.  The  lift  thjckness  deter¬ 
mines  the  computational  space  interval  (x)  which  is  one- 
quarter  of  a  lifl;  thiclness  for  a  given  computation. 

d. ,  Period  of  time  between  lifts.  The  number  of  days  (in  whole 

days)  between  lift'i  must  be  given.  The  number  may  be  varied 
to  allow  for  "irregular"  placement  scheduling,  or  may  be 
constant.  In  either  case  the  schediile  should  represent  in¬ 
sofar  as  possible  a  realistic  placing  schedule.  (The  time 
interval  T  used  by  the  computer  is  the  quarter-day  and  is 
not  affected  by  the  time  interval  between  lifts;  but  the 
n’jmber  of  quarter-days  M  is  four  times  the  number  of  days 
between  lifts.  M  should  be  furnished  computer  personnel.) 

t.  Adiabatic  temperature  rise.  Adiabatic  temperature  rise  data 
A0  are  obtained  in  advance  and  preferably  from  laboratory 
tests  of  the  concrete  mixture  to  be  used  for  the  specific 
job,  but  may  be  calculated  from  a  knowledge  of  heat  of  hy¬ 
dration  of  the  cement,  and  quantities  and  thermal  properties 
of  the  constituents  of  the  concrete.  A  separate  value  is 
needed  for  each  quarter-day.  The  computer  program  provides 
for  30  days  (120  quarter-days)  of  temperature-rise  data.  The 
increment  for  each  q\iarter-day  (A9)  rather  than  the  cumula¬ 
tive  rise  at  the  end  of  each  quarter-day  is  needed. 

f .  Thermal  diffusivity.  The  value  of  the  thermal  diffusivity  a 
is  determined  from  laboratory  tests  of  the  concrete,  and  is 
assumed  to  be  the  same  for  the  rock  foundation. 

£.  Foundation  temperature.  The  foundation-rock  temperature  (T^) 
must  be  specified,  but  may  be  assumed  to  be  equal  to  the 
air  temperature  of  the  first  ambient  period. 

5.  The  space  interval  X  ,  the  time  interval  T  ,  and  the  thermal 
dl ffusivlty  a  are  combined  into  the  single  dimensionless  constant  S 
for  use  by  the  computer.  Units  must  be  commensurate  in  the  relation  of 
these  variables,  given  by  the  formula  S  =  X  /(oT).  This  quantity  (S),  and 
not  its  components,  should  be  furnished  computer  personnel.  Since  T  has 
only  one  value  (0.25  day),  S  =  4x  /a.  lliough  not  likely  to  occur  so  in  the 
computation  for  a  dam,  S  must  not  be  less  than  unity.  The  A0  data  are 
given  in  a  manner  similar  ti)  the  following,  but  more  conveniently  in  tabu¬ 
lar  form:  1st  quarter-day,  6.2  F;  2d  quarter-day,  8.5  F;  3d  quarter-day, 
4.5  F;  4th  quartar-'^ay,  3.4  t?;...;  119th  quarter-day,  0.1  F;  120th 


quarter-day,  0.0  F.  Data  for  less  than  120  quarter-days  may  be  given, 
in  which  case  zeros  are  used  by  the  compuLcr  for  the  remaining  periods. 

6.  The  data  required  for  computer  so3.ution  of  a  problem  (see  para¬ 
graph  5  with  regard  to  S)  are,  therefore,  in  summary: 

a.  Ambient  (air)  temperature  (T  ),  to  hundredths  of  a  degree  F. 

b.  Concrete  temperatirre  (T^)»  to  hundredths  of  a  degree  F. 

c.  Lift  thickness  (4x),  feet  to  hundredths. 

d.  Period  of  time  between  placement  of  lifts  (M/4),  days. 

e.  Adiabatic  temperature  rise  (A0),  to  tenths. 

f.  Thermal  diffusivity  (a),  square  feet  per  day. 

g.  Foundation-rock  temperature  fT_),  to  hundredths  of  a  degree  F. 
Accuracy  to  hundredths  is  not  required,  though  all  data  may  be  submitted  to 
the  accuracy  shown  above.  Computations  within  the  computer  are  performed 
to  thousandths,  to  prevent  cumulative  error  which  might  occur-  in  isolated 
phases  of  the  computations,  but  the  computer  rounds  to  tenths  before  print¬ 
ing  computed  temperature. 


PART  II:  COMPUTER  PROGRAM 


7.  The  purpose  of  the  program  is  to  compute  temperatures  that  are 
expected  to  be  present  in  a  dam  during  construction  for  each  set  of  given 
conditions,  and  thus  to  provide  a  means  of  investigating,  prior  to  con¬ 
struction,  the  possible  effects  of  various  construction  practices  on 
temperature.  Temperatures  are  important  because  of  stresses  associated 
with  steep  temperature  gradients. 


The  Problem  and  Method  of  Solution 


8.  The  problem  was  to  write  a  program  for  solution  of  the  heat 
conduction  equation  for  one  dimension  modified  by  the  addition  of  a 
term  for  heat  generation  caused  by  chemical  action  of  the  cement  dtring 
hydration.  As  mentioned  earlier,  the  analysis  of  the  problem  and  method 
of  solution  were  described  by  Carlson,^  and  will  not  be  discussed 
further  herein.  Parameters  are  fully  discussed  in  Part  I  and  again 
with  illustrations  of  their  use  in  I’art  III,  where  the  accuracy  of  the 
method  also  is  discussed. 


Description  of  Equipment 


9.  The  program  is  written  for  a  basic  IBM-65O  computer  and  utilizes 
a  standard  utility  (8-word)  IBM- 533  panel.  It  is  computer-optimized 
by  the  Symbolic  Optimum  Assembly  Program  (SOAP)  and  operates  in  fixed 
point. 


Input-Output  Data  Fomnat 

10.  '‘‘  ..e  format  for  the  input  data  and  of  the  ou^-^nt  (answer)  list¬ 
ing  is  detailed  below,  and  the  specific  card  form  layout  is  given  in  the 
8-word  card  fo-vn,  fig.  i. 


Input  Data  Format 


Setup  Card  (One  Card) 


A0  Table  (1-120  Card) 


Lift  Cards  (l-50  Card) 


Identification  Card 


Detail  TenTperature  Card 


Maximum  Teii5)erature  Card 
(One  card  punched  each 
time  the  maximum  tempera¬ 
ture  computed  in  the 
quarter-day  interval 
exceeds  the  previous 
maxlmiim . ) 


1  -  Indent  -  This  word  may  contain  any  numeric  informa¬ 
tion  for  identification. 

2  -  Tpj  -  Tengjerature  of  the  fomdation  rock  is  fiven. 

3  -  S  -  The  value  of  "S"  is  computed  from  X^/(ofl')  smd 
given  herein. 

4  -  This  word  gives  the  value  of  M  vhen  the  con¬ 
crete  placing  schedule  is  constant .  "M"  is  the 

number  of  quarter-day  intervals  between  lifts.  If 
the  placing  schedule  is  variable  place  zeros  in 
this  word. 

Words 

1.  7, 

and  8  -  A9  table  routine  (ssune  for  all  cards). 

Word  3  ~  Quarter-day. 

Word  4  -  Heat  rise  value  A0  at  this  time. 

Word  1  -  "This  lift  number"  may  vary  from  one  to  fifty,  and 
is  the  identification  for  each  lift  to  be  placed. 

Word  2  -  Placing  temperature  of  the  concrete  for  thi.3  lift. 

Word  3  -  Air  temperature  for  the  period  of  exposure  of 
this  lift. 

Word  4  -  "M"  value  when  the  placement  scheduJ.e  is  varying. 

Word  4  of  setup  card  must  be  zeros.  For  constant 
placement  leave  this  word  blank,  that  is,  when 
word  4  of  setup  card  has  "M"  value  in  it. 

Output  Data  Format 

Word  1  -  The  Identification  word  of  the  input  data  is  given 
herein . 

Word  8  -  This  card  “is  identified  for  sorting  with  a  zero  in 
column  80. 

Word  1  -  The  number  of  the  top  lift  in  ,iaco. 

Word  2  -  The  number  of  thic  lift  for  which  the  following 
temperatures  are  given. 

Word  3  -  Temperature  in  upper  quarter  of  lift  at  this  time. 

Word  4  -  Temperature  of  middle  of  lift  at  this  time. 

Word  5  -  Temperature  of  lower  quarter  of  lift  at  this  time. 

Word  6  -  Teuiperature  at  bottom  of  lift  at  this  time. 

Word  7  -  Time  of  placement  of  this  lift  in  days  (quarter- 
day  Increments). 

Word  8  -  This  card  is  identified  for  sorting  with  a  1  in 
column  80,  Columns  7i-7^>  have  a  c.ard  count. 

Word  1  -  "Tins  lift  number"  in  vdiich  maximum  teuiperature  oc¬ 
curred  during  the  quarter-day  interval. 

Word  2  -  Time,  in  days,  at  which  maximum  temy  .’sture 
occurred. 

Word  3  -  The  maximum  temperature  that  occui'red  at  this 
time  and  lift. 

Word  8  -  This  i  ' -  identified  for  sorting  with  a  4  in 
column  0. 


Word 

Word 

Word 

Word 


WORO  1  WORO  2  WORD  3  WORD  4  |  WORD  5  WORD  6  I  WORD  7 
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Operating  Instruc''oions 


11.  Operating  instructions  for  the  program  are  c'/tlined  below; 
Program:  THERMAL  GRADIENTS  IN  MASS  CONCRETE  STRUCTURES  -  IO36D 


IBM-6^0  Console 


Storage  entry 

70 

1951 

70 

1951 

Programmed 

Stop 

Control 

Run 

Display 

Program 

Register 

Overflow 

Stop 

Error 

Stop 

3000(+)  (normal) 

990n(-)  (alternate  punch) 


Card  Input  fao  IBM- 333 

Program  deck  (7  per  card) 

Transfer  card  to  0200 
Setup  data  cai’d  (l  card) 

A9  Table  (1-120  cards) 

Transfer  card  to  0201 
Lift  data  cards  (I-50  cards) 

Start  Procedure 

New  Problem  -  I^ac.  pi’ogram  deck 

Trouble  stop  -  Check  for  error  stop  number  in  instruction  address  of 
program  register  and  compare  this  number  against  list¬ 
ing  of  error  stops  for  cause  and  remedial  action 
(page  9)- 

IBM- 533  Control  Panel  -  Standard  Utility  (8-word) 

Tabulating  machine  control  panel  -  Standard  Utility  (8-word) 

Special  operating  instructions  -  Sort  output  cards  on  column  80j  list  0, 

1  and  4  pockets  separately. 


Error  stops 


Instruc- 
Error  tion 
No.  Address 


Cause 


Transfer  to  location  0000. 
Possible  only  with  improper 
program  modification  or  fail¬ 
ure  of  650  to  read  program 
card. 

M  is  not  defined.  Check 
setup  card  or  lift  card. 

Address  of  Ti^  exceeds  storage 
in  rock.  Heat  flows  into 
rock  to  a  depth  equal  50 
lifts  of  dam  height. 

TLN  exceeds  TOPL.  Possible 
only  with  machine  error  or 
improper  program  modifica¬ 
tion 

I  exceeds  M.'  Possible  only 
with  machine  error  or  im¬ 
proper  program  modification. 

First  lift  card  read  by  pro¬ 
gram  is  not  first  lift  card. 


Console  switches  improperly 
set  up  for  alternate  piinch 
routine . 

Lift  card  read  is  not  proper 
card  to  be  read. 


Ti  address  not  equal  to  top 
lift  address.  Possible  only 
with  machine  error  or  im¬ 
proper  program  modification. 

L  address  not  eqml  to  lowest 
lift  address.  Possible  only 
witli  machine  error  or  im¬ 
proper  program  modification. 


Action 


Operator  reload  program 
deck,  if  same  stop  per¬ 
sists  ddmp  drimi  and  return 
to  programmer  for  action 

Return  to  initiating  of¬ 
fice  for  correction. 

Dump  drum  and  return  to 
programmer  for  action. 


Dump  drum  and  return  to 
programmer  for  action. 


Dump  drum  and  return  to 
pi^ogrammer  for  action. 


Operator  check  sequence  of 
lift  cards.  Place  cards 
in  proper  sequence  and 
transfer  to  0259  for  Read. 

Operator  check  for  sign 
equal  minus  and  switches 
equal,  70  1951  990(n). 

Operator  check  sequence  of 
lift  cards.  Reset  loca¬ 
tion  0204  with  value  equal 
to  next  lift  card  to  be 
read  minus  one  and  trans¬ 
fer  to  0639  for  Read. 

Diomp  drum  and  return  to 
pi-ogrammer  for  action. 


Dump  drum  and  return  to 
programmer  for  action. 


Un.'ess  program  has  been  modified,  error  stops  3;  8,  and  9  are  not 

likely  to  occur. 
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Flow  Diagram 

12.  The  flow  diag-ram  for  the  program  is  shown  in  fig.  2.  Modifica¬ 
tions  of  the  program  to  compute  slightly  different  problems  from  those 
discussed  will  require  use  and  understanding  of  the  flow  diagram. 

Alternate  punch  routine 

13.  The  program  with  normal  storage  entry  switch  settings  70  1951 
3000+  will  punch  one  detailed  temperature  card  for  each  lift  at  each  time 
interval.  The  number  of  output  cards  is  then  approximately  equal  to  but 
slightly  greater  than  that  computed  from  the  formula: 

s„=f  (H.l) 

Where  N  and  M  are  the  maximum  nimiber  of  lifts  to  be  computed  and  the 
number  of  quarter-days  between  lifts,  respectively. 

14.  The  niTiber  of  output  cards  punched  may  be  selectively  reducec" 
by  using  the  altamate  punch  routine  with  storage  entry  svritch  settings  of 
70  1951  990(n)  (-).  The  program  will  then  punch  ell  maximum  temperature 
cards,  but  only  the  selected  detail  temperature  cards  for  each  quarter-day 
computation,  (n)  in  the  switch  settings  indicates  the  number  of  intervals 
(quarter-days)  for  which  detail  temperature  cards  are  to  be  punched;  all 
inter'.nals  between  these  will  not  be  punched.  An  example  of  this  is  given 
below: 

Placing  schedule  (constant)  M  =  36  quarter-days  between  lifts 
(9  days). 

Alternate  punch  routine  in  use  with  (n)  set  equal  to  8,  storage 
entry  switch  70  1951  9908{-) 

Answer  cards  punched  out 

1.  fiaximum  temperature  cards  are  given  for  all  quarter-day 
time  interval  computations  as  they  occur  even  though  the 
detail  temperature  cards  are  not  punched  for  that  interval. 

2.  Detail  temperature  cards  are  gi^en  for  the  1st,  9'bh,  17th, 
25th,  and  33d  quarter-days  or  the  1st,  3d,  fth,  7th,  and 
9th  days  between  each  lift  placement. 

Computing  time 

15-  The  computing  time  of  a  problem  is  a  function  of  the  number  of 
lifts  to  he  computed  and  the  number  of  days  or  quarter-days  between  lifts 
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and  is  also  affected  by  the  use  of  the  alternate  punch  routine  which  re¬ 
duces  the  computing  time  slightly.  From  the  nature  of  the  program  no 
simple  equation  for  calculating  computing  time  can  be  given;  however, 
running  time  will  normally  vary  between  2  and  4  hours.  Problems  2  and  3 
of  the  Greers  Ferry  Dam  Report^  req.uired  3*2  and  2.8  hoixrs,  respectively. 
Other  routines 

l6.  A  Trace  routine  (311  TR)  and  a  ®even-Per-Card-Punch  routine 
(315  P?)  are  included  in  the  program  for  use  in  program  modification, 
checking  trouble  stops,  and  computing  long  problems  In  parts.  Reference 
is  made  to  IBM  Library  File  Nos.  1.4.001  and  1.2.002,  respectively,  for 
instructions  in  the  use  of  these  routines. 


Fig. 


2 


Flow  diagram 


Lk 


PART  III:  EXAMPLES 

Check  Problem 


17.  A  computation  of  temperatures  in  concrete  placed,  in  layers,  with 
concrete  placement  temperatiire  equal  to  the  ambient  temperature,  was  made, 
first  for  a  given  nimiber  of  days  between  lifts  (3  days),  and  then  for  two 
alternate  placement  schedules.  The  number  of  days  between  placement  of 
lifts,  constant  for  each  computation,  was  the  only  placement  condition  that 
was  varied  between  computations.  The  problem  was  completely  defined  by  the 
following  information  (temperatures  were  shifted  to  zero  for  convenience). 

Ambient  temperature:  0  deg  F. 

Concrete  temperatu’ce :  0  deg  F. 

Foundation-rock  temperature:  0  deg  F. 

Lift  thickness:  5  T't  (hence  1.25-ft  space  intervals  in  the 

computation) . 

Time  between  placement  of  lifts:  3  ^lays,  5  days,  and  7  days. 

Thermal  diffusivity:  1.00  sq  ft  per  day. 

Adiabatic  temperature  rise: 


Quarter-Days  AO,  F 


1  6.2 

2  8.5 

3  4.5 

4  3.4 


•  * 

119  0.1 

120  0.0 


The  proper  entry  of  these  data  on  WES  Form  No.  1053  is  shown  in  fig.  3* 

18.  During  a  computation,  temperatiires  are  printed  for  quarter-day 
intervals  and  for  three  positions  in  a  lift  and  at  the  boundaries  of  lifts. 
The  listing  for  the  l68th  quarter-day  time  interval  of  the  calculation  for 
3-day  placement  is  shown  in  ts.ble  1.  This  time  interval  corresponds  to  42 
days  and  l4  lifts.  The  temperatures  are  shown  plotted  in  fig.  4  (page  I8), 
where  also  are  plotted  similar  data  for  5-d.ay  and  7-day  placement. 
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Table  1 

Computed  Temperatures 


TOPL 

No. 

^1 

T 

2 

A- 

^4 

Time 

14 

14 

17.2 

26.8 

30.3 

30.1 

42.00 

14 

13 

30.0 

31.8 

34.0 

35.5 

42.00 

14 

12 

36.4 

37.4 

38.3 

39-0 

42.00 

14 

11 

39.6 

40.4 

41.1 

41.7 

42.00 

14 

10 

42.3 

42.8 

43.4 

44.0 

42.00 

14 

9 

44.5 

45.0 

45.4 

45.5 

42.00 

14 

8 

45.5 

45.7 

45.8 

45.8 

42.00 

14 

7 

45.8 

45.8 

45.8 

45.8 

42.00 

14 

6 

45.8 

45.7 

45.7 

45.7 

42.00 

14 

5 

45.7 

45.6 

45.5 

45.3 

42.00 

14 

4 

45.1 

44.8 

44.4 

43.9 

42.00 

14 

3 

43.3 

42.4 

41.4 

40.2 

42.00 

14 

2 

38.7 

37.0 

35.1 

33.0 

42.00 

14 

1 

30.7 

28.2 

25.6 

23.0 

42.00 

14 

— 

20.4 

17.8 

15.3 

13.0 

42.00 

14 

IM 

10.9 

9.0 

7.3 

5.8 

42.00 

14 

2M 

4.5 

3.5 

2.6 

1.9 

42.00 

14 

3M 

1.4 

1.0 

0.7 

0.4 

42.00 

14 

4M 

0.3 

0.1 

0.1 

---- 

42.00 

Note:  TOPL  is  the  nimiber  of  the  most  recently  placed  lift  of  concrete. 
Ti  is  1/4-,  T2  is  I/2-,  T3  is  3/4-lift  thickness  from  the  top  of  a 
given  lift,  and  Tj^^  is  at  the  lower  interface.  Time  is  in  days. 
Decimals  (as  above)  are  not  actually  shown  in  a  computer  listing. 
The  lata  shown  above  correspond  to  those  of  the  3-day  curve  in 
fig.  4.  The  computer  program  output  provides  a  complete  tempera¬ 
ture  history,  illustrated  by  the  above  data  for  one  time  interval, 
and  a  separate  listing  of  maximum  temperatures.  For  a  complete 
interpretation  of  output  listings,  see  reference  5- 


THERMAL  GRADIENTS  IN  MASS  CONCRETE  (IO36D) 


DATA  CARD  FORM 


Thermo/  R 0^(1  arch  and  Fia/d  Problem 
Submitted  by  Durahilil'^  Sech'on  No.  5  29  July  1953 


Program  constants 

Lift  thickness  (ItX)  5. GO  ft  Space  interval  (x). 

Thermal  diffusivity  (a)  1-00  ft^/dav  l/4  Lift  thickness  =  ft 

S  =  k^/a  =  G- 25 000 


Project 


SETUP  CARD 
Identification 


T 

17  ^  20 

s 

25  30 

M 

39  <+0 

XX.XX 

XX.XXXX 

XX. 

0000 

62500 

12 

WES  Form  No. 
Revised  Nov  1959 


Fif.  3 


Type  Cement 
^9 

23  26  37 _ !iO 

XXX.X 

0001  006  2 
0002  0085 
0003  OOU5 


[oopij.  005lf- 

0005  (ah.) 


0006 


HEAT  RISE  TABLE 


Type,  U ,  ^  hags /cu.  y</ 


23  26  37  ^0  23  26  37  *+0  23  26  37 


£B  BB 

40  37  40  37  40  37  4o  37  4o 


Key  punch  operator 

Use  load  card  format  below  for  each  card  of  heat 


0100 


rise  table  above » 


Word  1 

Word  3 

Word  4 

Word  7 

+ 

60  1954  1957 

+ 

20  XXXX  8000 

+ 

00  OCOO  XXXX 

•r 

19  0?0C  191;8 

Words  2,  5^  and  6  iray  be  blank 


I  I  I  {  I  I  I  I  i  I 

THE  AMBIENT  CONDITION  WAS  SHIFTED  TO  ZERO  IN  THIS  PROBLEM  POR  CONVEN¬ 
IENCE  THE  CORVES  COMPARE  PLACEMENT  AT  S-DAY,  S-DAY,  AND  7-DAY  INTER-  . 
VALS  THE  ADIABATIC  TEMPERATURE  RISE  USED  IN  THE  COMPUTATIONS  IS  ALSO 
PLOTTED  THE  PROGRAM  SEARCHES  POR  MAXIMUM  TEMPERATURE  IN  ONLY  THE 
FOURTH  TEMPERATURE  POSITIONS  IN  THE  CONCRETE  AS  SHOWN.  N  SLIGHTLY 
HIGHER  TEMPERATURE  MAY  OCCASIONALLY  BE  DETECTED  IN  ADJACENT 
POSITIONS.  T - 


I  IPT  NUMBER 


Fig.  4.  Electronic  computation  of  temperature  rise  of  concrete 

placed  in  layers 


A  Construction  Problem 


19.  Another  reproduction  of  Form  1053 ^  fig-  shows  similar  data 
for  an  actual  problem  for  a  construction  project.^  The  ambient  tempera¬ 
tures  shown  were  taken  from  a  curve  that  was  based  on  data  obtained  from 
Weather  Bureau  records.  Tlie  concrete  placement  temperature  (T^)  was  50  F 
for  this  problem.  The  results  of  this  computation  are  shown  plotted  in 
fig.  6.  Since  the  maximum  temperature  obtained  was  only  4  F  above  the 
foundation  temperature^  even  though  placement  of  concrete  continued  at 
9-day  intervals  throughout  the  summer,  the  advantage  of  placement  at  re- 
duced  teriperature  (refrigerated  concrete)  is  well  illustrated.  This 
problem  is  one  of  I6  such  problems,  representing  I6  combin8.tions  of  am¬ 
bient  s'onditions,  construction  practices,  cement  type,  and  thermal  proper¬ 
ties  c-f  concrete  wnich  were  investigated  to  determine  the  influence  of 
ihese  combinations  of  variables  or  of  the  separate  variables  on  the 


THERMAL  GRADIENTS  IN  MASS  CONCRETE  (10360^ 


DATA  CARD  FORM 


Tharmal  Research  anJ  Fiald 
Submit by  Durahillfg  SdcHon 


Program  constants 


Problem  Date  ^ 

No.  /£  2  July  1959 


Lift  thickness  7.50  ft  Space  interval  (x). 

Thermal  diffusivity  (a)  ft^/dav  1/4  Lift  thickness  =  1.875  ft 

s  =  4x2/a  =  9.765G>2 


SETUP  CARD 


Project 


Identification 


M 

39  40 


Grddrj  Farrg  Dom  00i03<o00l 2 


xx.xx 


7450 


xx.xxxx 


9765  Q 


M 

39  40 


XX.XX 


01  5000 


02  5000 


03  5000 


o4  !  (die) 


xx.xx 


7450 


7720 


7660 


7960 


6030 


8030 


6100 


8080 


6000 


7860 


7700 


7400 


7050 


0770 


0420 


0050 


S(b50 


5200 


4950 


LIFT  DATA  CARES 


20127  30  39  4o 


4040 


4570 


4130 


4150 


40  90 


4070 


4090 


4150 


rp  T  M 

20  27  ^  30  39  40 


Key  punch  operator  fill 
unused  columns  of  words 
1-4  v'ith  zeros  for 
standard  8- word  format. 
12  punch  columns  10,  20, 
30,  and  40. 


WES  Form  No. 
Revised  Nov  1959 


%o 


HEAT  RISE  TABLE 


Type  Cement  75  ~  25  hknd  of  PorHond  and  F/i/  ash 


r - 

23  26 

37  40 

23  26 

37  40 

23  26 

A0 

37  l^O 

23  26 

A0 

37  40 

37  4o 

37  40 

37  40 

37  40 

XXX.X 

0001 

0031 

0021 

0002 

004l 

0000 

0061 

0000 

0081 

0001 

0101 

0000 

0002 

00  Z  9 

0022 

0002 

0042 

0001 

0062 

0001 

0082 

0000 

0102 

0000 

0003 

0027 

0023 

0002 

0043 

0000 

0063 

0000 

0083 

0000 

0103 

OOOi 

0004 

OOZli- 

0024 

oOOi 

0044 

OOOI 

0064 

0001 

0084 

0000 

0104 

0000 

0005 

00 2  i 

0025 

0001 

0045 

0000 

0065 

0000 

0085 

OOOI 

0105 

0000 

0006 

0016 

0026 

0001 

0046 

0001 

0066 

OOOi 

0086 

0000 

0106 

0000 

0007 

00  i  2 

0027 

000  f 

0047 

0000 

0067 

0000 

0087 

0000 

0107 

0000 

0008 

0012 

0028 

0001 

0048 

OOOI 

0068 

0001 

0088 

0000 

0108 

0000 

0009 

0010 

0029 

0001 

0049 

0000 

0069 

0000 

0089 

OOOi 

0109 

0000 

0010 

0009 

0030 

0000 

0050 

0001 

0070 

0001 

0090 

0000 

0110 

OOOI 

0011 

0006 

0011 

0001 

0051 

0000 

0071 

0000 

0091 

0000 

0111 

0000 

003^ 

OOOG 

0032 

0000 

0052 

OOOi 

0072 

0001 

0092 

0000 

0112 

0000 

0013 

0006 

0033 

0001 

0053 

0000 

0073 

0000 

0093 

OOOi 

0U3 

0000 

0014 

0006 

0034 

0000 

0054 

0001 

0074 

0000 

0094 

0000 

0114 

0000 

0015 

0005 

0035 

0001 

0000 

001^ 

0001 

0095 

0000 

0115 

0000 

0016 

000^ 

0036 

0001 

0056 

0001 

0076 

0000 

0096 

0000 

0116 

0000 

0017 

0003 

0037 

0000 

0057 

0000 

0077 

0000 

0097 

OOOi 

0117 

OOOi 

0018 

0003 

0038 

0001 

0058 

00(1 

0078 

0001 

0098 

0000 

0118 

0000 

0019 

0003 

0039 

0000 

0059 

0  00 

0079 

0000 

0099 

0000 

0119 

ooco 

0020 

0002 

0040 

OOOi 

0060 

0001 

0080 

0000 

0100 

0000 

0120 

— 

Key  punch  operator 

Use  load  card  format  below  for 


each  card  of  heat  rise  uable  above. 


Word  1 


60  1951+  1957 


Word  3 


20  XXXX  8000 


Word  4 


00  0000  XXXX 


Word  7 


19  0202  Y)'/o 


Word  8 


31  0002  i9;;3 


Words  ?,  5,  and  6  may  be  b.lank 


Fig.  5  (Continued) 
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Fi.g.  6.  A  piediction  for  a  construction  project  (Greers  Ferry  Dam, 

problem  12) 

predicted  thermal  history  of  the  concrete  in  Greers  Ferry  Dam.^ 


Accuracy 

20.  The  accuracy  Indicated  for  each  item  of  data  in  paragraph  6  is 
satisfactory  for  a  computation.  However,  before  preparations  for  the 
Greers  Ferry  computations  (preceding  paragraph)  proceeded  very  far,  it  be¬ 
came  apparent  that  the  accuracy  of  certain  data  for  repetitious  use  by  the 
computer  would  have  to  be  increased  by  several  decimal  places  beyond  the 
accuracy  that  had  been  indicated  in  brief  manual  computations  and  also  in 
electronic  computations  for  the  "clieck  problem"  (paragraph  17)»  The  value 
of  S  must  be  to  five  significant  figures  (or  six  if  S  S  lO.OOOO),  as 
shown  in  completed  Foims  1053  (fig-  5)  and  1013  (fig-  l)>  and  other  data 
should  be  as  shown  in  these  forms.  Although  it  was  impracticable  to  check 
a  complete  problem  by  manual  computation,  a  satisfactory  method  of  checking 
the  general  accuracy  of  results  of  a  computation  was  devised  and  used.  T>ie 
procedure  was  as  follows;  a  finite  increment  equation,  based  on  the  orig¬ 
inal  differential  equation, "^’3  vas  applied  as  shown  in  fig.  7.  In  this  way 
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close  approximation  checks  could  he  made  at  any  questionable  point  on  a  tem- 
peratxire  curve.  Several  discrepancies  were  isolated  by  this  method,  and 
necessary  adjustments  to  tlie  program  were  made.  It  is  believed  that 
sources  of  sy&bematic  error  were  reduced  to  a  practical  minimum. 


0.10 


^U6.25* 

_ z 

1 

■—213. 

1 

i 

25* 

AO 


L\V  I  '*0.15  I 

—  =  (1.4.))  - 

[  (7.5)2  J 


♦age  of  OAM  (DAYS 
AFTER  placement 
OF  FIRST  LIFT). 


=  0.00766  P  PGR  DAY,  OR 
0.43  F  IN  56  DAYS  AT 
CHECK  POINT  0  IN 
AGREEMENT  WITH  THE 
DIFFERENCE  OETWEEN 
COMPUTED  CURVES. 


=  0.00384  F  PER  DAY, 
OR  0.257  F  DROP  IN 
67  DAYS  AT  ChEICK 
POINT  0  IN 
AGREEMENT  WITH 
THE  DIFFERENCE 
BETWEEN  COMPUTED 
CURVES. 


SINCE  MANUAL  COM* 
PUTATION  OVER  A 
LONG  PERIOD  OP 
Time  was  IMfRACTI- 
CAL,  T  ^E  METHOD 
ILLUSTRATED  IN 
THIS  FIGURE  WAS 
USED  T  CHECK  THE 
ACCURACY  OF  THE 
COMPUTED  TEMPERA* 
TURE  CHANGES. 

_ j _ 


"  TBMPERATURB  DURING  EXPOSURE  OF  LIFT,  IN  DEGREES  F 
I  69.7  I  73.0  I  76.1  \  77,8  |  79.0  {  80.0  |  30.7  |  8J.0  |  81.0 


3  4 

LIFT  NUMBER 


Fig.  7*  Check  of  accuracy  of  computed  temperature  changes 


Concluding  Remarks 

21.  The  cost  of  electronic  computation  of  temperatures  in  a  dam  is 
reasonably  estimated  to  be  less  than  one-tenth  the  cost  of  manu/il  computa¬ 
tion.  The  cost  of  the  electronic  computation  of  temperatures  at  Greers 
Ferry  Dam  vas  moderate  compared  vith  the  total  cost  of  uhe  customary 
thermal  investigations  that  are  made  for  most  suen  projects.  Ma>iy  of  the 
expensive  thermal  tests^  such  as  those  for  thermal  diffusivity  and  adia¬ 
batic  heat  rise^  have  been  developed  principally  because  the  results  of 
these  tests  are  used  in  temperature  gradient  computations;  however^  because 
of  the  _"ost,  manual  temperature  gradient  computations  have  not  been 
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undertaken  very  extensively  in  the  ^st.  These  therms,!  properties  tests, 

t 

which  hitherto  have  had  only  limited  usefulness,  provide  the  .data  tt^at  are 

needed  for  a  complete  thermal  analysis,  as  is  'uell  illustrated  hy  the  re- 

5  ' 

suits  of  the  computatj.ons  for  Greers  Ferty  Dam.  It  appears  at  the  present 
time  that  no  serious  emulative  error  is  involved  in  the  computer  program, 
and  that  for  comparative  purposes  temperatures  predicted  hy  means  of  the 
eJ.ectronic  computations  are  entirely  valid.  It  is  recommended  that  ‘at 
least  several  basic  computations  be  made  for  all  major  projects  in  the 

I 

future  * 

) 
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